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A 0(h ) Cubic Spline Collocation Method for Quasilinear 

Parabolic Equations * 



1. Introduction . Consider the quasilinear parabolic equation 

(1.1a) c(x,t,u) u - u = f(x,t,u,u ) , 0 < x < 1 , 0<t<T 

t XX X — 

(1.1b) u(0 , t) = b (t) , u(l,t) = b.(t) , 0 < t < T , 

• O 1 

(1.1c) u(x,0) = g(x) , 0 < x < 1 



Let A = {0 = x q < x^ < • • • < = 1} be a partition of I = [0,1] , with 

I. = [x. - x.], h. = x. - x. .. , and h = max (h.: 1 < i < N} . Then define 
l l-l, i l l l-l l — — 



and 



n k (J) = V:{V is a polynomial of degree < k on J } 



n = (V: V £ IL (I.) , 1 < i < N> . 

k,A k i — — 



For -!<_-£<_ k-2 let 

s(A,k,£) = n k>A n c l (i) 

be the space of piecewise polynomials of degree < k (order = k) on A with 

continuity t . Note that S(A,k,£) has dimension d[S(A,k,£)] = kN - (£+ 1) (N-l) . 

In this paper we shall be primarily concerned with S, = S(A>4,2) , the usual 

n 

cubic spline space on A . 



This research was supported in part by the National Aeronautics and 
Space Administration (NASA) . 
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Recently Douglas and Dupont [ 10, 11 ] have studied collocation procedures 

for (1.1) based on the spaces S(A,k,l) for k > 4 . Their main result is 

that collocation at the images of the k Gauss-Legendre points in each sub- 

interval 1^ yields uniform errors of order 0(h) and superconvergence 

2k — 4 

results at the knots {x^} of order 0(h ) if the solution of (1.1) 

k+2 

u e H (I) . These estimates (but not the analysis) are essentially the 
same as those of deBoor and Swartz [ 7 ] for ordinary differential equations. 

The analysis of [11 ] is based on certain discrete innerproducts as is the 
analysis presented in this paper. 

Several authors [4 ,16 ,20 ] have studied collocation techniques for 

ordinary differential equations using smoother spaces S(A,k,£) with Z > 2. 

k-2 

The general result obtained is that the convergence rate is 0(h ) , a 

suboptimal rate of convergence for such spaces. However, the procedures of 

Russell and Shampine [20 ] will provide 0(h ) convergence (and superconvergence 

2k- 6 

at the knots of order 0(h )) for k > 6 if the collocation takes place at 

the images of the k-2 Lobatto points on each subinterval. Hence, it is ex- 
pected that these procedures can be extended to parabolic problems through 
careful mimicing of the arguments in [ 11]. These procedures will be studied 
in a later paper. 

2 

In [ 8 , 19] cubic spline methods with 0(h ) accuracy have been studied 
for linear versions of (1.1). Also, in [18 ] a cubic spline collocation pro- 
cedure for the heat equation has been proposed (but not analyzed); for a partic- 

4 2 

ular explicit time discretization 0(h + (At) ) convergence obtains. However, 

this procedure is essentially the standard explicit finite difference method 
for the heat equation and the high accuracy does not readily generalize to more 
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difficult problems or other time discretizations. 

In [ 2, 3] a variant of the usual cubic spline collocation method yielding 

4 . 

0(h) convergence rates for nonlinear ordinary differential and quasilinear 
parabolic equations was studied. In this paper we describe the high-order pro- 
cedure and provide continuous time estimates for (1.1). Two types of boundary 

4 

collocation will be considered, yielding uniform estimates of 0(h ) and 
7/2 

0(h ) respectively. In a subsequent paper we shall investigate the effect 

of various boundary collocation techniques for high-order smooth spline approx- 
imations to (1.1). 

It should be noted that the particular approximation used here is essentially 
the same as that of Daniel and Swartz [ 9 ] for two point boundary value problems. 
The procedures were developed independently, the derivation of [2 ] preceeding 
that of [ 9 ]. The finite-difference method discussed by Hirsh [13 ] can 
be interpreted as a cubic spline method, and as such, it is quite similar to the 
present technique. 

This paper has four parts. In § 2 the basic notation of the paper is developed. 
Some discrete innerproducts for cubic splines are then defined and studied. The 
basic approximation technique used here is developed in § 3 from consideration 
of a simple two-point boundary value problem. In § 4 the main results are pre- 
sented (Theorem 4.1). 
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2. Notation . We shall use the standard notation [ 14] for L^(I) spaces 

m 2 

and Sobolev spaces H (I). In particular, | |v| | „ = (v,v) with 

IT (I) 



(f 



,g) = / f(x)g(x)dx . Also, H 1 (I) = {V e H 1 (I) : v(0) = v(l) = 0} , 

J I ° 

. Also, we use W m (I) = (v: D^v is abs. cont. 



with v 



hJ (I) 



Dv 



l 2 (i) 



m oo 

0 2 j < m-1 , D v e L (I)} with 



v 



= l || A | 



W (I) l < m L~(I) 

The spaces L^[0,T; X] are defined as usual for normed linear spaces X [ 14 ] . 

OU 

If v £ L (I) is defined on A , then write v_^ = v(x^) and v = (v q , , • • • , V ^ T ) 



Let 



v = max v. 

1 1 1 n OO 1 I 1 

l 0<i<N 



Define the difference operators Vv. = h.^(v. - v. _) , Vv. = Vv, , and 

r ill l-l l l+l 

2 — 

A = VVv^ . In case v = v(x,y) , denote the differences with respect to a 



particular variable as 



(V^v) (x^y) = h i 1 (v(x i ,y) - v(x i _ 1> y)) ; 



etc . 



In the following A is uniform; i.e., A = {x^ =ih:0^i<N}. Then 

define the discrete innerproduct 
h N 

[v,w] - 2 £<Vl”i-l + V i W i > 

l=l 



with norm | v 1 2 = [v,v] 2 . Also, let < v,w > = h \ v i W i anc * l v l = < v,v >2 * 

N 1=1 

Additionally, let < v,w ] = h £ v.w. 

i=l 1 1 



Recall the summation by parts formula: 



(2.1) < Vv,w > = - < v,Vw ] + v n w n - v^w q 
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The following results are easily established for cubic splines. Let 

S° = S, n H 1 (I) . 

h ho 



Lemma 2.1 If v,w e S, , then 
h 

2 4 2 

(2.2a) - < v'',w > = (v',w') + Y 2 (v'',w'') + (v" ' ,w" ' ) - B(v' ' ,w' ) 

2 4 2 

(2.2b) - < v" + j 2 A v” ’ w > = (v'.w’) + Iso ' I 2 B ( v ”’ w ’> » 

h 4 u 2 

(2.2c) [v’.w’J = (v\w’) - ^20 (v" , ,w”') + ^2 {B ( V ”’ W ’> + B(w",v')} , 

.2 

(2. 2d) [v",w"] = (v",w") + f- (v"',w'") , 

where 

B(v,w) * v X 7 w x7 - v w 
N N o o 

Note that v,w e (I) is not necessary for (2. 2d). 

o 

PROOF : We prove (2.2a-b); the remaining results are similar. Recall the 

corrected trapezoidal rule 



d 2 

y <Kx) dx = [<Kc) + <f> (d) ] — [(f>' 



d . (d-c) 5 , (4) r . n 

J + —^q — <P (S) > 5 e (c,d) 

C 



Applying this rule one interval at a time and summing yields 

2 __ 4 2 

(v' ' ,w) = [v” ,w] - ^ < Vv” , Vw ] + Tm (V " ,v' ” ) - ^ B(v' ’ ,w' ) 



180 



12 



for all v,w £ . Summation by parts and v,w e (I) imply that 

2 4 2 

(2.3) (v",w) = [v”,w] + ~ < A 2 v",w > + Yso (v’ ",w"') - yy B(v",w') 



which is (2.2b). It is easy to show that 

(2.4) (v' ' ,w") = B(v' ' ,w' ) + < A 2 v' ' ,w > ; 

hence, (2.2a) follows from (2.3) and (2.4). 
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To apply these results, we need the inverse relations (not assumptions 
for splines) . 



Lemma 2.2 


t 21 ] 


If V 


€ n k A Wlth 


h/ min hj_ < a , then 








K. , i± 


l<i<N 


(2.5) 


1 |dv| 


1 a 


< Ch _1 | | v | | 


~ > 






L q (I) 


1 _ 1 


L q (I) 


(2.6) 


i|v|| 


L q (I) 


£ Ch q P || 


v | | , l<p<q<«> 

L P (I) 




Let 


|v | 

1 1 8 


= max | v (x) | 
x=0, 1 


. Then for v e S° , 
h 


(2.7) 


|b(v’ 


',v’)| 


1 2 |v"l 3 |v’ 


l 3 < Ch" 2 I |v|| 2 1 

d H X (I) 

o 



Hence, by (2.2c) 

(2.8) | v | 2 + |v' | 2 < C| | v | | 2 

H (I) 
o 



Since B(v ft ,v T ) is not definite, the left sides of (2.2a) and (2.2b) 

(with v = w) are not norms equivalent to the H^(I) norm on S° in general. 

o h 

Of course, if v is periodic, B(v Tl t v ! ) = 0 and the forms in question are 

actually equivalent to the H^(I) norm. 

o 

It is not generally true that |v| and | |v| | ? are equivalent; 

L Z (I) 

however, it is the case that 



(2.9) I v l C | | v | | 2 

L (I) 

It is true that 

(2.10) C, ||v|| 2 £ | v | 2 + h 5 |v" | 2 £ C.J | v | | 2 

1 l/(I) L (I) 
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For this note first of all that the exponent of h is correct by Lemma 2.2. 



Also, if 


i ,2 , 5, .. ,2 


= 0 , 




l v l + h |v | 3 


then 




o 

II 

•H 

> 


0 < i < N 


and 




O 

II 

- *H 
> 


i = 0,N 


It is then clear that 


v = 0 . Hence, 


lllvlll H 


on S° and (2.10) 

h 


follows . 





' 2 + h 5 |v 






is a norm 



We shall also use the following notation 



, tr 2 

< v,w > ^ < v + — A v , w > 






for mesh functions v and w . 
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3. A Two-Point Boundary Value Problem . In this section we consider a cubic 
spline approximation to the two-point boundary value problem 



(3.1a) 

(3.1b) 



u"(x) = f(x) 



x e I 



u(0) = b , u ( 1 ) = b 
o 1 



,N 



It is well-known [ 4 ] that collocation at the knots {x_^K_q 

; i.e., finding e such that 



in 



(3.2a) 

(3.2b) 



U f 1 (x.) = f (x.) 

Cl 1 



0 < i < N 



U (0) = , U (1) = b 1 

c o c 1 



has a convergence rate 0(h) (and no better) in general. However, defining 
U 6 S h by 



(3.3a) 

(3.3b) 



U f 1 (x. ) = f (x. ) - — f * 1 (x. ) , 0 < i < N 

l l 12 i 



U(0) = b , U(l) = b- 



leads to the following results. 



Theorem 3.1 Suppose u e wV) is the solution of (3.1) and U e is 
defined by (3.3). Then the following estimates hold for e = u - U : 



(3.4a) | |D j e| | < Ch 4_j | |u| | , 0 < j < 3 . 

L (I) W (I) 

The following super convergence results are also valid for 1 < i < N : 

(3.4b) | e .' | < Ch 4 | |u | | , 

1 W°(I) 

(3.4c) |e' | £ Ch 4 ||u|| 

1 2 w b (i) 

00 

(3. 4d) | e' '(£..) | < Ch 3 | |D 5 u| | ; j =1 , 2 , 

LJ L (I.) 

(3.4e) £Ch 2 ||D 5 u|| ^ 

1 2 L (I.) 



9 



3-2 



where 
(3. 4f ) 
and 
(3. 4g) 



‘i-k = ( Vl + X i )/2 



5.. = x ^ + (-1) J , j = 1,2 

1J 1 * 2YT 

,6+k , 



Additionally, if u e W (I) for 0 < k < 2 , then 

(aV')« = U'. ' + o(h 4+iC | |D b+ u| 

L W (I) 



(3.4h) 



U* ! + -^r (aV ') . = U'.' + o/h 4+k | |D 6+ u| | ), 

i 12 i i \ M L“(I)/ 



for 1 < i < N-l 



Proof : Expand e' ' about x. . on I. to obtain for 0 ^ T ^ h : 



(3.5) 



(x + T) .fid. + 

^ i-1 } \ 12 2 2 / i-1 V 6 12 / i-1 

+ °(» 4 l|f lv ll . ) • 

\ L (I,)/ 



( V' 



Then it is straightforward that 
(3.6) 



J' e"(x)p 2 (x) dx | <_ Ch 5 ||u| 



I. 

l 



w 6 (i ± ) 



1 £ i £ N , 



for any ^2 € ^2 A bounded independently of h • 

Let G^(x;0 be the Green’s function for v T ' = g on I subject to 

v(0) = v(l) = 0 , and define G (x ;£) = ( “ G J(x;5) . Since 

1 y 3x o / 

G (x; •) e W 1 (I) and G.(x.; •) e II for 0 < i < N , j = 0,1 , 

O J 1 Z. y Q 



we have 



|D J e(x )| = I fc Ax. ; Oe"(0 d? | 

I J 1 

N f 

1 l \J G. (x ; C)e"(0 dC | 

n=l I J 



(3.7) 



N 



1 Ch l | | u | 



n=l W (I ) 

n 



< Ch u 



W 6 (I) 
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The stability results of [ 22 J and (3.7) imply 



(3.8) 



I D J e I 



' 00 
L (I) 



Ch 



u 



1 6 

W (I) 



Integration of (3.5) from x. to x. l yields (3.4c). Estimates 

1 “ 1 

(3.4d-e) follow immediately from (3.5). 

Since for 1 < i < N-l and 0 < k < 2 

2 

(3.9) (A 2 U"). = [A 2 (f - ^2 f”)] i = C' + o(h 2+k | |D 4+k f | | m ) 

L (I ) 

estimate (3.4h) is established, and the proof is complete. 



We now consider defining 




W e by (3.4h) neglecting 

More precisely, define W e S, 

h 



(3.10a) 


W"(x) = f (x) - Y 2 f"00 




h 2 7 


(3.10b) 


W! • + 77 (AW"). = f. 

1 12 11 


(3.10c) 


W(0) = b , W(l) = b- 

o 1 



the 

by 



The following results then obtain . 

Corollary 3.2 Let u , U be as in Theorem 3.1. Define W e by 



(3.10). 


Let z = U-W 


and e = z + e 


= u - W . Then 


(3.11a) 


IIdLi 


1 . 1 Ch 4+k 1 


|d 6 *II. 






L (I) 


L (I) 


(3.11b) 


lh'" 


II. 1 ch 3+k | 


|l> 6 +k u|l„ 






L (I) 


L (I) 


if 


U € W 6 +k (I) 


, 0 < k < 2 . 





Furthermore, all the inequalities of (3.4) hold with e replacing e 
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and the norm of u on the right side of each inequality being changed 



to 



w 6 (i) 



Proof : Equations (3.10a-b) yield an (N+l x N+l) linear system for W. 



(3.12) 

where 



A = 



12 



A W” = f = f - tx (fl' > 0 , * * * , 0 , f ' ' ) 



0 

0 



12 ' o 



12 0 0 0 



1 10 1 
0 1 10 



0 0 



• • # 



• • • 



N 



• • • 



0 0 



0 

0 

0 



10 1 
0 12 



™*1 3 

Since A is diagonally dominant with j | A | | , W f T (hence W) 

oo C. ^ 

is uniquely defined by (3.11). By (3.4h) U f f satisfies 



(3.14) 

where 

and 



A U f T = f + 6 , 



6 = (0,6 1 <S N _ 1 >0) 



| 6. | < Ch 4+k ||u 6+k | 



L (I) 



, 1 < i < N-l 



Subtracting (3.12) from (3.14) yields 
(3.15) A z" = 6 
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from which 



(3.16) 



£ fl l« 



< Ch 4+k ||D 6+k ul 



L (I) 



Since z is piecewise linear 

(3.17) III" || _ £ I I I Z" I 

1 ( 1 ) 



< Ch 4+k ||D 6+k ul 



L (I) 



Finally, the Green's function representation leads to for j = 0,1 



(3.18) |D : ’z(x)| = | f G (x ; Oz" (O d£ | 1 | | G (x ; •) 

I J J 



w ou 

L (I) L (I) 



which along with (3.17) establishes (3.11a). The remainder of the Corollary 
follows from homogeneity in h and the triangle inequality applied to 
e = z 4- e 
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4_. Continuous Time Estimates . In this section we consider the continuous- 
in-time approximation to the solution of (1.1) by cubic spline methods. 

Define U : [0,T] by 

2 

(4 ’ la) (“xx + 12 c(U >"xxt ) < x ’ t > ‘ (V U)0 t - V U)U xt - F h (U 'V) < x >*>> X ' O’ 1 

2 

(4.1b) (c(U)U - [U + jt A 2 U ] ) (x. ,t) = (f (U,U ) ) (x. ,t) , 1 < i < N-l , 

\ L XX X XX /I \ X / 1 

(4.1c) U(0,t) = b (t), U(l, t) = b_(t) , 

o 1 



(4. Id) U(x,0) - g(x) = "small" , x e I . 



In (4.1a) we have used (supressing (x,t)) 



h 2 2 



(4.2a) A^*) = c(*) - YJ D h £c(<f») J 



.2 1 

(4.2b) B h (*) = \ [c(*)] 



(4.2c) 



h 2 2 

f(<M x ) - 12 D h I f (4>. 4> x ) J 



where 



(4. 2d) D^[<H(x) 



1 (V x «1 



2 (7 X W 1 



2 (, x *h 



2 (V x *Vl • 



X = 0 



x = 1 , 
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and 






| *) x - 


(A 2 <M 2 , x 


= 0 


(4.2e) 


l>](x) 




2<4 X 


(A x ^^N-2 * 


x = 1 . 


Note 


that for 


Ip £ W 2+k (I) , 


0 < k < 2 




(4. 3a) 


nt 

/-N 


D i 1) 

X 


- D l [« - 


°(h k || Dk+ %|| 


L (I) / 


(4.3b) 


< <♦> - 


D 2 ^ 

X 


- D l ['] - 


o(h k ||D k+2 ^|| 


. )• 

L (I) ' 


In 


the following we shall assume that c and f are smooth functions of 


their arguments with c 


subject to 


the bounds 




(4.4) 


0 < m _< c 


(x, t,<j>) <_ H < °° 








x e I , t 


£ [0,T] and <$> 


e R = (- 00 , 00 ) 





The choice of this particular approximation is motivated by the results 



of §3, specifically Corollary 3.2. Later U(x,0) will be chosen to provide 
the desired 0(h ) convergence rates. It is also possible to use a different 
collocation procedure at the boundary, namely 



(4 . la) ! 



(c(U)U -U )(x,t) = f(x,t,U,U ) , x = 0,1 ; 

\ t XX / X 

7/2, 



however, the analysis here will provide only 0(h ) rates of convergence. 

The analysis of (4.1) will proceed along the same lines as that in [ 11 ] 
and will employ the discrete inner products of §2 . Before beginning the 
error analysis, we establish the existence and uniqueness of the solution of 
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(4.3) 



(4.1). For this, we consider the equivalent matrix formulation based 

B-spline basis (V- , V 0 , . . . , V, T , _ } on the knot set 

12 N+3 

I 0 = x. (1< i < 4), t 4+ . = X. (1 < i < N-l), t n+3+ . = 1 , (1 < i 1 

N+3 



see [ 5 , 


6 ]. Let U(x,t) = £ 

J = 


a. (t) V . (x) 
1 3 3 


(4. 5a) 


^(a)a'(t) -Ja - J 


(a) , t e 


(4.5b) 


a i(t) = bg(t) 


’ V3 (t) ' 


(4.5c) 


II 

/-N 

O 


given , 


where 


a(t) = (a 1 (t),a 2 (t). 





c( l \ V i, ( x ±) ) ( x ±) » 1 < i < N-l 

k K K J 



(4.6a) = 



l 



12 c( l “k V x i )>V j' <x i> - A h ( l o k V k ( Xj 
k J k 



+ B h ( ^ a k V k (x i ))V i (x i ) ’ 1 = °’ N 

1 r J 



(4.6b) 1^]^ 



v l (x i> + To (a2 v -')( x i) » 1 ^ 1 s N- 

J -L^ x J 



-r. • (x.) 



, i = 0,N 



(4.6c) 






f ( I * l \ v ^ x i)) » 1 S i < N- 

k k 



_F h ( ^ “k \ (x i> ’ l \ = 0,N 

k k 



on the 



4) }; 



))v j (x ± ) 



1 



1 
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The assumption that f is Lipschitz continuous with respect to its 
last two arguments implies that ^(a) is likewise Lipschitz continuous. 
Thus, the local existence (in time) of the solution to (4.1) will be es- 
tablished in case 



(4.7) 


✓'“N 

P 

II 

o 


and 


B 1 P N+3 ° 








N+3 


implies that 


6=0. Let 


4> (x) 


= £ a, V (x) and suppose that 








k=l k k 



N+3 

iKx) = £ 3, V (x) with 8 satisfying (4.7). Then 

k=l K ^ 



(4.8a) <Kx-£) =0 , 0 < i < N , 

and 

(4.8b) {c(<j>)ip + 2D^ [c(<f>)]^ }(x.) =0 , i = 0,N . 

xx n x 1 

By the standard cubic spline identities [15 (4.8) is equivalent to 



(4. 9a) 



(4.9b) 



(4.9c) 



i P (x. .) + 4i|) (x i ) + ip (x. .. ) = 0 
XX x— 1 XX 1 XX 1+1 



{c(<t0 - Y dJ [c (4>) ] } (0) , l J xx (o) - 3 



{c(<)>) + ~ D* [c(<|>)]}(l)^ xx (l) + ~ 



, 1 < i < N-l , 

[ c (4>) ] = 0 » 

[ c (4>) ](l)ij/ (x ) = 0 . 

XX N— 1 



Assuming that 

(A. 10) |~ [c(x, t, <f>) ] | £ L for x £ I , t £ (0, T], <J> £ R , 

we find from (4.3a) that {D^ [ c (4>) ] } (x) is bounded for x = 0,1 . Thus, 

n 

for sufficiently small h , (4.9) corresponds to a diagonally dominant, homo- 
geneous linear system for ^ (x-^) • Hence, = 0 and 6=0. 

XX XX 
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Lemma 4.1 


If (4.10) holds then for h sufficiently small there exists a 


unique U £ 


solving (4.1) for t € (0,T] . 


Lemma 4.2 


if 


(4.10a) ' 
and 


1 [ C (x> t J <j>) ] | £ L , X £ I 9 t £ (0 9 T ] 9 (j) £ R 9 


(4.10b)' 


||^ [f(x,t,<M)]| £ L , x € I , t € (0 , T] 9 <M € R, 


then for h 


sufficiently small, there exists a unique solution 



of (4.1) with (4.1a) f replacing (4.1a) for t e (0,T] . 

Proof : Similar to the above. Just differentiate (4.1a) ? with respect to t 

to obtain an analogue of (4.9) which is diagonally dominant for h small enough 

We now turn to the convergence analysis of (4.1). Note that (4.1b) is 
equivalent to the discrete Galerkin formulation 

(4.11) < c (U) U ,V > - < U ,V > = < f(U,U ),V > , V € S u . 

t xx A x h 

For the analysis define the comparison function W : [0, T] -* by 



(4.12a) 

(4.12b) 


h 2 

W x x ( *" t) * u x X <X ’ t) - T2 u xxxx (,tlt) 1 ’“°' 1 > 

(W xx + l2 A x V 1 "! 1 ' 1 ■ u x* <X l- t ) ■ 1 S 1 < N-l , 


(4.12c) 


w(0,t) = b 0 (t) , W(l,t) = b 1 (t) 



4 

Note that Corollary 3.2 implies that W (W fc ) is a 0(h ) approximation to 
u (u t ) . 
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Let z = W - U £ , e = u — W , and e=z+e=u-U . Our 

plan is to estimate z in terms of e and then to bound e using the bounds 

on z and the triangle inequality. In the following analysis, we shall often 

2 2 

require the inequality ab < ea + (l/4e)b for a,b < 0 , any e > 0 . 

From (4.12b) and (1.1) we find 

< c(W) W , V >-< W , V > = - < c(W) e ,V > 

t xx A t 



(4.13) 



+ < [ c (W) - c(u)]u ,V > + < f(u,u ),V > , V e S, 

t x h 



Subtract (4.11) from (4.13) and apply the assumed smoothness of c and 
f to obtain 



< c(U)z t ,V > 



< z ,V > = < [ c (W) - c (U) ]W , V > - < c(W)e ,V > 

XX A t t 



+ < [c(W) - c(u) ]u t> V > 



+ < f(u,u ) - f(U,U ), V > 



(4.14) 



= < c u zW ,V > - < c(W)e t ,V > 



- < c e u , V > . 
u t 

+ < f e+f e,V> 

u u x 

x 

* * 

+ <f z+f z ,V> , 

u u x 

x 



where the partial derivatives c ,f ,f are evaluated as required by the mean 

u u u 

x 

value theorem. Now use Cauchy-Schwarz , the boundedness of the derivatives of 
c and f , and the trivial inequality mentioned above with V = z^_ : 



(4.15) < c(U)z t ,z t > - < z xx »z t > A - C ( | z | 2 + |z x | 2 + | Te | ) + 6 | z | 2 , 6 > 0 
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Here 

(4.16) 



I Te | 2 = |e| 2 + |ej 2 + |ej 2 



If 

(4.17) 



u,u fc e L 2 [0,T;W 6 (I)] , 



Corollary 3.2 implies that 

(4.17b) / | Te | 2 dt < Ch 8 (| | u | | 2 „ , + 

0 ' L Z [0, T;W° (I) ] 



C L 2 [0,T;W 6 (I)] 



Choose <5 so that n = m-6 > 0 and use (2.8) to obtain 



(4.18) 



n|zj 2 ~ < z xx ,z t ^ - C( l | 2 1 + l Te l 2 ) 

t xx t A H 1 (I) 



To complete the estimate, it is necessary to consider the boundary terms 
(4.1a). A straight-forward computation using (4.12a) and (4.3) yields 



{ “ xx + 12 c(W)W X *t Kx ’ l:) ' {A h (u)u t ‘ B h (u)u x t' V u>u x )}<x ' t) 

, 2 



(4.19) 



where 



12 {c(u)e xxt }(x ’ t} 



R^(x>t) , x — 0,1 , 



h 2 '2 



(4.20a) (x, t) = ~ J2 ^ " 2R ^ ^ [f (u, u^) ] } (x, t) 



Note that 



if I h = [0, x 3 ] u [x N _ 3 ,l] , 



(4. 20b) 
and 

(4.20c) 
then by (4. 3) 
(4. 20d) 



c (x, t ) = c(x,t,u) e W (1^) 



f(x,t) = f(x,t,u,u x ) £ W (I h ) , t e (0, TJ 



£ Ch 4 (||D 4 c ||^ <v + ll#ll_. \ W 



( V 
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In the sequel, we assume that 
(4. 20e) = 



D x C M - co + ll D * f M 

L [0,T;L (I h )] x 



< OO 



L [0,T;L (1^)] 



Subtract (4.1a) from (4.19) to obtain 
. 2 



z + — c(u)z 
xx 12 xxt 



xt 



(V u > - V u) h 

- (B h (u) - B h (U))u 

- (B h (u> (z xc + ; xt >) 

(4.21) - (F (u, u ) - F (u, U ) ) 

n x n. x 

- < F h (u "V - F h (U ’V ) 

■ 12 c(u) \ x t + \ 

We now estimate the terms on the right side of (4.21). The treatment of all 
but one of the terms is somewhat rough. 

IV") - \ (u) l 3 i l A h ( "» - v w >l 8 + l*h (w > - \ <u >i 3 



(4.22a) 



2 

= Y2 ( [c(u) " c(W)] l 3 + l D h [c(w) 

■ T5 0 D h l D h [ < zl 0 

* CCII |e|| | + || M|| ) 



:(U) | 



Similarly, 




(4.22b) 


| B h (u) - B h (U)| a < Ch(| | |ej | | ra + || |zj | |J 


(4.22c) 


' F h (“"Vls- c(ll|; J | l- + ll|z x 


and 




(4. 22d) 


| p h ( u »U x ) - F h (U,U x ) | 3 < C(|| |e| 11^+ | | |z| | | e 
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Recalling (4.10) and (4.3a) with k = 0 , 



Thus , 
(4. 22e) 



lV u >l s i Ch ‘ 



l B h (l,Kz x t + Vis i Ch VJa + 



Now multiply (4.21) by z > integrate in t , and apply (4.22) 

XX L 



to obtain 



1 I z 

2 xx 1 3 



r t 

2 2 ft. 2 -21 

a (t) + yh / |z xxt l 9 dt < ch J 



| z| I I 2 + I I |z III 2 + h 4 | z I 2 

i i ioo i • i x 1 1 1 xt 1 3 



(4.23) 



+ l Te l 3 + 



i\i lb 



+ 2 l z Ja <0) • 



where y > 0 



(4.24a) | Te | 2 = | | |e| | | 2 + 



, , , I 2 x u^lr I 2 j. I~ I 2 

e + h I e . \ r, + e . L. 

x 1 1 ,QO \ 1 xt 1 3 1 xxt 1 3 



(4.24-1,) / | Te | 2 dx <_ Ch 8 J (||u|| 2 6 +||u|| 2 6 ] dx „ 

0 3 0 ' W (I) C VT(I) / 



if C 4 . 17 a) holds. 



It is clear from (4.21) and the bounds (4.22) that 

2 



lz Ja- 1" Vth + cl 11121 



+ |||z III + h 2 | z I 

x 00 1 xt 1 8 



(4.25) 



l T ®U + 15,1a) 



, t £ (0 , T ] . 



Use of (4.23) and (4.25) permits the completion of the estimate (4.18). 



From (2.2b) 



(4.26) 



- < z z > = — 

xx J t A 2 dt 

« 2 



z|I hJ(I> + llo H z x X kI^ 2 (i) 



- — B(z ,z ) . 
12 xx xt 
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Observing (via Lemma 2.2) that for any e > 0 



12 i b <s*-vi 



(4.27) 



< e | |z | | 2 + Ch|z | 2 , 

t L z (I) xx 3 



and adding (4.27) to both sides of (4.18) yields 



I |2 1 d / I I II 2 

n Ki + 2 



IKJI 2 2 

(!) X L (I) 



(4.28) 



i c ( I l z l | 2 X + h l 

\ hJ«) 



z | 2 + | Te | 2 
xx 1 3 1 1 



+ e z 



l 2 (D 



Integrate (4.28) with respect to t , apply (4.25) to bound the h z 
term, and apply Lemma 2.2 to produce 

„t 



XX 3 



2 dt +-l||z || 2 (t) £ c/jlUII 2 , 

'0 L 2 H L (I) 01 n 

o o 



n/ U t l 



+ Te 



(4.29) 



+ h ( | Te | 2 + |rJ 2 )| dt 



c t , z 

+ Kh J~ I z dx 

0 1 xxt 1 3 



2 \ £ 



+ ( E + K h ) /IU t 



dx 



L (I) 



+ C | | z | | (0) . 

H (I) 
o 
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Now add (n 
of (4.23) . 




(4.30) 



+ K)h 5 




z 



XX t 




to (4.29) using the estimate 



Thus, proceeding as in (4.29) 



dT + 



ilMl\ <t) < c/jiuil 2 , + 

H (I) 0 ( H (I) 



Te 



+ h ( | Te | g + iRjg) dt 



+ 



(e + K*h 




2 

l 2 (i) 



dx 



+ c I l z l | 2 1 (0) 

H A (I) 

o 
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By (2.10) and Gronwall's inequality with e and h sufficiently small, 
we obtain our basic estimate 



(4.31) 



\H 2 2 2 + IUII 2 

l/[0,T;l/(I)j 



L [0,T;H o (I)] 



1 C } | 
n ( 



2 

3 



Te | + h 1 1 Te | ^ + 

2 






dx 



+ c|| Z |r (o) 

H (I) 
o 

If (1.1) is smooth enough that (4.17a) and (4.20e) hold, then the use of 
(4.17b), (4.20d), and (4.24b) in (4.31) implies the estimate 



(4.32) 



z t l | 2 2 2 + INI 2 

L [0,T;L (I)] 



L [0,T; H (I)] 
o 



< Ch 






L [0,T; W (I) ] 



L 2 [0,T;W 6 (I) ] 



<) 



+ C| | z| | (0) 

K (I) 

o 



There are a variety of interpolation schemes which produce U(x,0) such 
that | | z|| I (0) = 0(h^ f ) . See [15] for several possible choices. Here we 

H o (I) 

choose the interpolation studied in [9]; namely, define U(x,0) by 



(4.33a) 

(4.33b) 



U(x i ,0) = g(x i ) 



0 < i < N 



U (x,0) = g"(x) “ g (iV) (x) , x - 0,1. 

xx 12 



It is then easy to see (compare with Corollary 3.2) that if g e W (I) 



(4.34) 



xx 



L (I) 



(0) Ji Ch ||g|| 6 



W (I) 
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Our final estimate is then 



(4.35) 



zJI 2 2 + IW 

* L [0, T; L (I) ] 



< Ch 



L [0,T;H o (I)] 



V ' 2 



L [0, T;W (I)] 



+ N*N.6_ +K h) 



W (I) 



after using the fact that for <f> e L [0,T;X] , <j>(0) e X , 



< C 



L [ 0 , T; X] 



C L 2 [0,T;X] 



( 0 ) 



The further assumption that u e L [0,T;W^(I)] , and Corollary 3.2 imply that 



(4.36) 



oo oo 



< Ch u 



L [0,T; L (I)] 



L [0, T;W (I)] 



Using (4.35), (4.36), the triangle inequality and the embedding of H X (I) in 

o 

< C I 1 4 > I I .. , <(> e H^(I) , we obtain the uniform 

L° (I) ~ Hhl) 

. o 

estimate 



OO 

L (I) J i. e. , 



< Ch u 



L [0,T; L (I)] 



L [0,T;W (I)] 



(4.37) 



* t l I 2 



L [0, T; W (I)] 



+ ^) 



Theorem 4.1 Suppose that c, c , f, f , f are uniformly bounded mde- 
u u u 

x 



pendently of their arguments and that (4.10) and (4.20b,c) hold. Then for 
h sufficiently small there exists unique U solving (4.1) and (4.33). If 
u , the solution of (1.1), satisfies 



(4.38) 
then 

(4.39) 



ue L°°[0,T;W 6 (I) ] , u t e L 2 [0, T;W 6 (I) ] , 



lu-ui | w £ Ch ( Mul 

L [0, T; L (I)] 



+ U 9 , 

t L 2 [0,T;W 6 (I)] 



L [0, T; W (I)] 

+ K h)- 
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We now show that the use of the simpler boundary collocation (4.1a) 1 

7 / 2 N # 2 °° 

yields 0(h ) estimates; suboptimal in L (I) or L (I) norms, but 

optimal in H^(I) • This loss of accuracy is not believed to be actual; 

rather, it is just a function of the particular analysis employed. We 

proceed as before with U defined by (4.1), (4.1a) 1 and W by (4.12). 

Then the error analysis is unchanged through (4.18). Note that for x = 0,1 



(4.41) 



c (W) W - W = f(u,u ) + e 

t xx x xx 



where 

h 2 

(«.«> ««(*.<> ■ T 2 u xxxx <x ' t) ' x * °- 1 

Subtract (4.1a)* from (4.41) and use the boundary values to find 



(4.43) 



-z = f(u,u ) - f(U,U ) + e 
xx x x xx 



= f (z + e)+e 
u x x xx 
x 



Thus , 



(4.44) 



ifc r**j ^ 

B (z ,z ) = B (f z ,z ) + B (f e + e , z ) 
xx xt u x’ xt u x xx xt 

x x 



= — - B(f z ,— z ) + B(f e + e , z ) 
2 u x dt x u x xx xt 

x x 



Integrate by parts 




+ 



B 




e + 
x 



xx 



z 

X 



) 



t 

o 
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Now assume that 



(4.46) 



3 2 f 
3t3u 



x 



(x, t,(f>,ijj) I £ L < °° , xel.te [0,T],<J>,4i e R. 



Then 



B(Z xx’ Z xt )dT l - (c + l Z x l 3 dT + Ch ! I T ^3 dT 



(4.47) 



+ (C + e h 1 ) |z x | 2 (t) + Ch | Re | 2 (t) 



where, 



+ C(h 1 | z x lg + h | Re | 2 ) (0) , 



(4.48a) | Te | 



l\>3 + ^ 6 xx ^ 3 + l e xtl 3 + ^xxJs 



and 



(4.48b) Re 



l : J» + l ; J? 



Integrating (4.18) with respect to t and applying (4.26) and (4.46) yields 



/ |zj 2 dT + ^ 



h\i) 

o 



(t) £ c / (| | z | | 2 1 



+ Te 



H o a) 

+ h 3 1 Te | 2 ) dx 

2 r 2 

+ (Ch + eh) Jjzjg dT 
+ (C h + e h) | z J ^ ( t) 



x 1 3 



+ Ch 3 | Re | 2 (t) 



+ C 



(h|z x l 9 + h I Re | 3 



1 >> 



h o (D 
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Apply Lemma 2.2 several times and take h so small that 




* 

E ) > 0 . 



Then 





I Te | 2 + h 3 1 Te | 2 ) dx 



(4.50) 



+ 




max 

0<T<t 



I Re I ^ )+ IUII 2 . (0) 

3 H (I) 



Gronwall's inequality implies that 



f T \ z I 

•'n't 1 



2 dr + | | z | | 2 



L [O.TJH (I)] 



(4.51) 



<C f ( 



I Te | 2 + h 3 |Te| 2 )dt 



+ C 



(h 3 ( max | R e| 2 )+ | | z | ] 2 (0) ) 

\ o<t<T w m f 



H o (I) 



Note that 



n cox . 3 | ~ . 2 _ _h7 , , 2 3,~ .2 _ h 7 , .2 

' e xx^3 144 ' U xxxx^3 * ^ e xxt^3 144 ^ U xxxxt ‘ 3 



Hence* under assumption (4.17a) and the choice of U(x,0) given in (4.33), we 



obtain 

/l* t | 2 * + ||.|| 2 „ ! 

0 L [0,T;H o (T)] 



(4.53) 



Ch 



u 



t 



2 

L 2 [0,T;W 6 (I)] 



+ 



2 

w 6 (i) 



+ 




U 

xxxxt 



2 

a 



dx + max 
°<r<T 



u 

xxxx 
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This estimate leads easily to the optimal 

1 7/2 

H q (I) norm; however, it implies only 0(h ) 

00 

or L (I) norm. Of course if 

(4.54) u (x,t) = 0 , x = 0,1 , t € [0,T] , 

xxxx 

4 

the estimate becomes 0(h ) . Similarly, if u is periodic, i.e., 

(4.55) u(0,t) = u(l, t) , 0 £ j £ 2 , t e [O.T] , 

4 

then making U similarly periodic results in 0(h ) estimates since 

B(z ,z ) = 0. Note that in this case (4.1a) is replaced by (4.55) as 
xx xt 

applied to U . 



0(h ) estimate in the 

2 

estimate in the L (I) 



Theorem 4.2 Under the assumptions of Theorem 4.1 , with (4.10) f and 
(4.46) replacing (4.10) and (4.20), let U be defined by (4.1) with (4.1a) 1 
replacing (4.1a). Then, 



(4. 56a) 
(4.56b) 



I I u-U 



II 1 < Ch 3 ( 


Mull r + 


U 1 


1 t 00 - \ 


00 r 6 / 


M t 



L [0,T;W (I)]> 



u-U 



Il»[0,T;l"(I)] - Ch ,l|u| 



+ u 



L [0,T;W (I) ] 



fcl ' L 2 [0 ,T;W 6 (I) ] 



+ CcJ ^ i max | u 
' 0<t<T 



xxxx 1 3 






If, in addition, (4.54) holds 
(4. 56c) | I u-U I I 

il l I oo oo — 

L [0 ,T;L (I)] 

Furthermore, if (4.55) holds, and 
(4.56c) obtains. 



Ch 



.6 



t 1 1 2 6 

L [0 ,T;W° (I) ]/ 



L [0,T; W (I)] 

U is also required to be periodic, then 
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